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Abstract. In thin (A < few A;) collisionless current sheets in 
a space plasma like the magnetospheric tail or magnetopause 
current layer, magnetic fields can grow from thermal fluctua- 
tion level by the action of the non-magnetic Weibel instabil- 
ity (Weibel T959J. The instability is driven by the counter- 
streaming electron inflow from the 'ion diffusion' (ion in- 
ertial Hall) region into the inner current (electron inertial) 
region from where the ambient magnetic fields are excluded 
when released by the inflowing electrons which become non- 
magnetic on scales smaller than the electron gyroradius and 
< few A e . It is shown that under magnetospheric tail con- 
ditions it takes ~ 20-40 e-folding times (~ 10-20 s) for the 
Weibel field to reach observable amplitudes |b\y| ~ 1 n T. In 
counter-streaming inflows these fields are predominantly of 
guide field type. This is of interest in magnetic guide field 
reconnection. Guide fields are known to possibly provid- 
ing the conditions required for the onset of bursty reconnec- 
tion ([Drake eTaTj |2006l [Pritchett| |2005l [20061 |Cassak et al.| 
|2007| >. In non-symmetric inflows the Weibel field might it- 
self evolve a component normal to the current sheet which 
could also contribute to reconnection onset. 

Keywords. Weibel fields in thin current sheets, Weibel 
thermal level, Guide field generation, Magnetospheric sub- 
storms, 



1 Introduction 

In this communication we investigate the self-consistent gen- 
eration of a so-called magnetic guide field in a thin collision- 
less current layer that initially lacks the presence of any guide 
field. For this we render responsible the Weibel instability 
HWeibel[[T959) . 



Guide fields are believed - and have been shown by nu 
merical simulations (see e.g. |Drake et al. 2006 Pritchett 
[2005] [2006] [Cassak et al.||2007| and others) - to be of prime 



importance in collisionless reconnection. Their presence in a 
thin collisionless current sheet is not evident. In a collision- 
less plasma the state of magnetisation is determined by the 
history of the plasma. External magnetic fields are unable 
to penetrate it; they have to be present in the plasma from 
the very beginning. It is thus of vital interest in reconnection 
to understand whether guide fields can arise in collisionless 
current sheets. 



Space observations in situ ( 


Fujimoto et al. 


1997; 


Nagai et 


|aL]|1998||0ieroset et al.||2001 


, Runov et al. 2003 ; Nakamura 



lations (Ramos et al. 



2002; Drake et al. , 2003] |2005| l unam- 



biguously prove that magnetic reconnection in a collisionless 
space plasma proceeds under the following two conditions: 

- that the current sheet that separates the oppositely di- 
rected ('anti -parallel') magnetic fields ±B to both sides 
of the current becomes 'thin enough'. Under 'thin 
enough' it is understood that the effective half-widths 
| A < Xi of the sheets fall (approximately) below the 
ion inertial scale-length A, = c/uji (with c velocity of 
light, oj{ =e(iV/eomj)2 ion plasma frequency, e el- 
ementary charge, N plasma number density, mj ion 
mass). The real limit on the width is not precisely 
known. Observations (Naka mura et al.[ |2002[ |Vaivads| 
et al. 2004| |Nakamura et al.| |2006[ |Baumjohann et al] 



2007) suggest that it can reach values A < 4Aj. This 



condition is commonly (for example in Retino et al. 
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2007 1 Sundqvist et al. 2007) taken as sole indication 
that one is dealing with collisionless reconnection. 



- that plasma flows in into the sheet from both sides along 
the normal to the current at a finite (though low) velocity 
±V n , explicitly considered only in treatments of 'driven 
reconnection' . Numerical simulations usually take care 



2 



R.A. Treumann, R. Nakamura, and W.Baumjohann: Weibel fields in thin current layers 



of this condition by starting the simulation with a pre- 
scribed reconnection configuration either assuming an 
initially present X-point or by locally imposing a tem- 
porary parallel electric field or finite resistance for suf- 
ficiently long time in order to ignite reconnection. 

The first condition implies that ions in the sheet become 
inertia-dominated and thus are demagnetised, while electrons 
remain magnetised. Since magnetised electrons are tied to 
the magnetic field, the magnetic field in the ion inertial re- 
gion is carried along by the electrons causing Hall current 
flow (as was realised first by Sonnerup 7979}. 

From the second condition it is clear that reconnection 
cannot proceed continuously on time scales shorter than the 
inflow time r- m ~ A/V n ~ iewXi/V n . If reconnection turns 
out to be faster, it will necessarily be non-stationary and 
probably pulsed. 

These just represent necessary conditions still being in- 
sufficient to describe the onset of reconnection. A collision- 
less mechanism is missing so far that either demonstrates, 
in which way the electrons become scattered away from the 
magnetic field in order for letting the oppositely directed 
magnetic field components slide from the electrons and re- 
connect, or that forces reconnection to occur in some other 
way. Prime attention in this respect is attributed to the ac- 
tion of guide fields, i.e., a magnetic field component in the 
direction of current flow. 

Magnetic guide fields have two implications which in a 
thin current sheet with unmagnetised ions affect basically 
only the electrons: 

- that the centre of the current sheet is not free of mag- 
netic fields, as would be the case in the Harris current 
sheet model, and 

- that the sheet current J possesses a field-aligned com- 
ponent Jii which, when exceeding a certain threshold, 
undergoes instability and may cause either anomalous 
dissipation via generation of anomalous resistance or 
viscosity, or produces localised wave fields like solitons 
and (electron) holes. Both resistance and guide-field- 
aligned electric fields violate the electron frozen-in con- 
dition and thus contribute to (localised) reconnection. 

Quite generally, any guide-field-free current sheet model 
that scatters electrons away from the magnetic field turns out 
to be in trouble because of the following reason: 

Assume that the electrons have transported the magnetic 
field from both sides some distance across the ion inertial 
region into the current sheet. During this transport the mag- 
netic field lines have become bent locally thereby enhancing 
the magnetic tension. However, magnetic fields are massless 
and therefore have no inertia. Once the electrons enter the 
electron inertial region in the centre of the sheet, they de- 
magnetise and release the magnetic field lines. When this 
happens, the field lines react to the magnetic stress and will, 



hence, rebound to their initial straight magnetic field con- 
figuration, stopping moving further in for meeting their op- 
positely directed partner-field lines and merging with them. 
This behaviour will keep the internal section of the current 
sheet free of magnetic fields. (The stationary kinetic Har- 
ris equilibrium does not cover this effect; it just prescribes 
a tanh-shape for the magnetic field and does not distinguish 
between ion and electron inertial regions.) 

Clearly, the rebounding field lines will again become 
loaded with electrons when snapping back into the 'ion diffu- 
sion' Hall-current region being unable to completely regain 
their initial positions. Thus, magnetic field lines accumulate 
near the boundary of the electron inertial region until - for 
continuous further plasma inflow (e.g. in driven reconnec- 
tion scenarios) - the accumulating external pressure becomes 
strong enough to push the newly created magnetic wall fur- 
ther in. 

In this respect it is of interest that very few of the magnetic 
profiles across a presumably reconnecting current sheet at the 
Earth's magnetopause (cf., e.g., Haala nd et al.| |20"04| |2005| 
Paschmann et al. 2005) or also in the magnetospheric plasma 
sheet ( |Runov et all ]2003 ) exhibit the canonical Harris tanh- 



shapes; most of them show the evolution of either a mag- 
netic plateau across the current sheet centre or regions along 
the profile where the magnetic field spatially undulates. The 
latter behaviour is conventionally attributed to waves pass- 
ing along the boundary (Paschmann et al. |2005 i referring 
to related changes in the boundary normal. This seems rea- 
sonable; part of it might, however, also be explained quite 
naturally by the above effect, in particular when there is tan- 
gential flow along the current layer - as is the case at the 
magnetopause. 

Reconnection takes place only, when the oppositely di- 
rected magnetic field lines are brought into direct physical 
contact to (partially) cancel each other, i.e. coming so close 
that the field-lines touch each otherQThis, however, requires 
some mechanism that acts even deep inside the electron in- 
ertial region on scales < fcwA e = c/w e ~ few km (< 10 km 
in Earth's tail plasma sheet and even ~ 1 km in the mag- 
netopause region; with u; e = e-J 'N / 'eom e w 56.41-v/iV krad/s 
angular electron plasma frequency, m e electron mass, and N 
is in cm -3 ) transporting the magnetic field further in. It is 
not known in which way non-magnetic electrons could cause 
further inward transport of the magnetic field. In quantum 
solid state physics such a transport is related to the cele- 
brated integer and fractional quantum Halls effects (see, e.g., 
Ezawa 2002 for an almost complete account). However, 



due to their high temperatures T 3> Tp (with Tp Fermi tem- 



*A precise formulation of this condition requires the definition 
of a field line cross section which cannot be given in classical 
plasma physics but requires a quantum electro-dynamical approach 
(R.A.Treumann, The quantum picture of reconnection, in prepara- 
tion, 2010). 
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Plasma inflow 




Plasma inflow 



Fig. 1. Sketch of a homogeneous thin collisionless current sheet with plasma inflow from both sides and central electron inertial region. 
The magnetic field piles up at the boundary of the electron inertial region giving rise to a broad (of size of a few electron inertial lengths 
A e = c/uj e ) field-free central current layer of sheet current J which is crossed by two (symmetric counter-streaming) electron flows which 
are Weibel unstable. On the right the pile-up deformed Harris profile is shown. It exhibits two shoulders in the regions where the magnetic 
field piles up at the boundaries of the broad field-free central electron inertial region, where the Weibel instability may generate magnetic 
guide fields. 



perature) space plasmas are classical, and - presumably - no 
similar possibility opens up here. 

Thus the question arises whether a collisionless current 
sheet can by itself cause a guide field to evolve in its centre 
which circumvents the application of an external guide field, 
re-magnetises the demagnetised electrons in the centre of the 
current sheet and thus allows transport of the lobe magnetic 
field into the centre of the current sheet, contact of oppositely 
directed fields, and finally makes onset of reconnection pos- 
sible. 

In this Letter we demonstrate that the Weibel instability 



(Weibel 1959 1 is capable of doing this sufficiently fast. The 
newly generated Weibel guide-field could become compara- 
bly strong, a non-negligible fraction of the undisturbed ex- 
ternal field. As expected it is directed along the current 
layer and thus indeed plays the role of a guide field, suggest- 
ing that under certain conditions collisionless narrow cur- 
rent sheets self-consistently generate weak guide fields of 
strength B g /Bo < 1 where Bo is the magnetic field strength 
in the inflow region. The current sheet then becomes locally 
capable of producing conditions under which reconnection 
takes place. 



2 Weibel scenario 



The Weibel instability (Weibel 1959| l is known to produce 
stationary magnetic fields under conditions when the plasma 
exhibits certain anisotropies in flow and/or temperature. It is 
thus of interest to investigate the conditions for its excitation 
in a collisionless thin current sheet. This instability is driven 
either by electrons or ions with the ion instability being much 
weaker than the electron instability. 



The original proposal by Weibel ( 1959) referred to a tem- 
perature anisotropy in the unmagnetised electron distribution 
providing the free energy for a stationary (very low frequency 
uj ~ 0) magnetic instability. Various variants of this instabil- 



ity have been investigated in th e past two decades (|Yoon and 
|Davidson| [19871 |Yoonl [2007al |Medvedev and Loeb[ |19991 
ICalifano at ai~] 2002; Sil va et al.| [20021 |Fonseca et al.||2003] 



Achterberg and Wiersma 2007 1 both in the non-relativistic 



and relativistic domains, unmagnetised and magnetised, cold 
and hot plasmas, and for temperature as well as beam in- 
stabilities, with the beams effectively faking a temperature 
anisotropy. Application of this instability was mostly in- 
tended in either laser (inertial) plasma fusion or the violent 
conditions present in astrophysical systems. As we will show 
below, reconnection provides one of the simplest and most 
interesting classical applications of the Weibel instability. 

Figure [T] sketches the collisionless reconnection site, con- 
ventionally called the 'ion-diffusion region' (even though 
there is no diffusion). The ions just become non-magnetic 
here. The frozen-in electrons have sufficient momentum to 
continue their inward E x B-drift motion transporting the 
magnetic field to the centre of the thin current sheet. Close to 
the centre the electrons become demagnetised and release the 
magnetic field which snaps back as explained above. In the 
narrow electron inertial region of size of few A e perpendic- 
ular to the current sheet the collisionless inert electrons still 
maintain their inward velocity V = ±14 z (with 14 < V n ) on 
both sides of the current sheet. Since the electrons are com- 
pletely collisionless, the two flows pass across each other 
without (direct) interaction thereby realising a non-magnetic 
counter-streaming electron beam configuration which ac- 



cording to Weibel (1959) may become electromagnetically 
unstable. 
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The Weibel instability generates a non-oscillating (u> ~ 0) 
transverse magnetic field b with k-vector about perpendic- 
ular to the electron beams (k± ^> fcii, subscripts refer to the 
direction of V). In space plasma this instability is manifestly 
non-relativistic, and as usual one might suspect that it is too 
weak to cause any susceptible effect. We will, however, show 
below that its effect is not negligible. Investigating instability 
we work in the fluid approximation of cold (Tf, -C T e ) beams 
of density Nb, also for simplicity assuming that the plasma 
is cold as well, even though the centre of the reconnecting 
current layer contains a denser N > Nb thermal plasma of 
temperature T — T e +Ti (to which we will return when esti- 
mating the thermal fluctuation level of the instability). Under 
unmagnetised plasma conditions in the presence of beams 
the Weibel instability grows. So far it has been investigated 
mostly in the relativistic limit. In reconnection in space plas- 
mas inclusion of relativistic effects is unnecessary. The in- 
flow Mach numbers are small, and beam velocities Vf, -C c 
are small as well. 

The two symmetric counterstreaming electron beams in 
the central electron-inertial region have roughly same den- 
sity Nb, a condition that applies primarily to the magneto- 
spheric tail-current sheet; at the magnetopause the beams are 
non-symmetric. At very low frequencies the electromagnetic 
dispersion relation factorises (see also the more elaborate pa- 
pers by |Yoon and Davidson||1987{|Achterberg and Wiersma| 
2007 and many others where non-symmetric, thermal and 



kinetic corrections are given and also external magnetic fields 
are applied, with the expected result that the latter act stabil- 
ising on the evolution of the instability). 

In slab geometry applicable to the geomagnetic tail cur- 
rent sheet, the factor describing plane electromagnetic fluc- 
tuations of frequency cj s» becomes 



D^D 



XX " ZZ 



ID, 







(1) 



where Dy are the components of the dispersion tensor 
D(w,k) = k 2 c 2 \ — e(a>,k) (with plasma dielectric tensor 
function e(w,k)). Under the assumed symmetric conditions 
D xz = 0, and the dispersion relation simplifies to 
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Fig. 2. Conditions inside a local electron inertial region. The ex- 
ternal magnetic field accumulates at the boundary of the inertial 
region. Electron flow is perpendicular to the magnetic field. In- 
side the inertial region the electron components form two oblique 
beams. The Weibel instability generates a magnetic vortex of field 
bw with wavenumber k _L Vj, perpendicular to the beam config- 
uration. These vortices close in the current plane (x,y) and have 
wavelengths of the order of ~ A e b. The sense of the magnetic vor- 
tices alternates along the current direction. Shown is the part of 
the lowest order vortex where the field is along the sheet current J. 
At a phase of tt/2 along the current the Weibel fields point along 
the ambient magnetic field causing wobbling of the current layer. 
Under non-symmetric beam conditions the wave vector and Weibel 
fields may become tilted against the current sheet. This could give 
rise to the local production of a normal magnetic field component 
b z which may possibly initiate reconnection. 



where n is the refraction index, k wave number, u> wave 
frequency, and s \ij the susceptibility tensor of species s 
the only surviving component of which, in a symmetric 
electron/electron-beam plasma configuration, is given by 



Xz 



1- 



(3) 



The subscripts e and b indicate background and beam pa- 
rameters, respectively, ujb is the beam plasma frequency for 
symmetric beam density N± = Nb/2, and in the last term 
the (negligibly small) neutralising background-ion contri- 
bution to the plasma frequency is taken into account for 



correctness in the electron-to-ion mass ratio term m e /mj. 
When the background plasma is at rest, the 'wave' becomes 
non-oscillating with to = ±i^fw (otherwise when the plasma 
moves at velocity V the wave frequency will be Doppler 
shifted by the amount k • Vrj, a case that may be realised un- 
der magnetopause conditions, there with Vn being the mag- 
netosheath flow velocity tangential to the magnetopause). 

Solving the above dispersion relation for 7w > yields 
the growth rate of the non-evanescent growing mode 




2N 



rrii 



(4) 
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Again, N is the quasi-neutral background density. Maximum 
growth rates 



7W,r 



(5) 



are obtained when the second term in the braced expression 
becomes small, which is the case at relatively short wave- 
lengths 



fcA efc »[l + (2A/A fc )]5>l 



(6) 



The Weibel instability will thus lead to comparably small- 
scale magnetic structures populating the beam-electron in- 
ertial range of size A ~ 27rA e fc = 2irc/ui{,. This condition 
of maximum growth seems to restrict the range of unstable 
wavelengths to very short scales. However, since the beam 
density is substantially less than the density of the ambient 
plasma in the centre of the current sheet, one has A e b > A e , 
the condition that k < A e is thus not in contradiction with the 
condition of maximum growth. 

The unstably generated magnetic structures are not purely 
transverse, however, because kc/u > 1. They are of mixed 
polarisation, containing a longitudinal electric 'wave' field 
component eg | k which in the symmetric beam case is small, 
however. The main transverse electric field component ex is 
along the electron beams. Since the 'wave' magnetic field 
b satisfies the solenoidal condition b • k = 0, k _L Vf,, and 
wb = kxeT, the Weibel magnetic field lies in the plane that 
is perpendicular to the plasma inflow into the current sheet 
(cf. Figure |2j. 

In an extended two-dimensional current sheet the Weibel 
magnetic field is directed either parallel or anti-parallel to 
the sheet current J = Jy while being confined to the electron 
inertial zone around the non-magnetic centre of the current 
sheet. Such a field has the required properties of a guide 
field. Since it cannot be stronger than the magnetic field B 
in the external inflow region, the Weibel instability gener- 
ates weak guide fields satisfying bw/Bo < 1. The above 
smallness condition on its wavelength implies that the Weibel 
guide field forms a comparably short-scale wavy magnetic 
structure along the sheet current, thereby structuring the cur- 
rent sheet magnetically in the direction of the current flow. 
In addition it has a transverse electric field eT||V& compo- 
nent along the electron beam inflow direction that is confined 
to the electron inertial region, while the longitudinal electric 
field is in the direction et _L Vfc. 

The case that the inflow is not homogeneous either along 
the current or transverse to it is of particular interest. Then 
the wave vector kj^ rotates in the (x,y)-plane perpendicular 
to the inflowing electron beam. The condition that the mag- 
netic field be free of divergence forces the Weibel magnetic 
field to form closed magnetic vortices in this plane. Such 
a magnetic vortex causes a magnetic asymmetry along the 
sheet current by periodically amplifying and weakening the 
external magnetic field on one side of (above or below) the 
current layer. This produces a spatial oscillation around the 



symmetry plane of the current layer of wavelength of the 
Weibel magnetic vortices. 

3 Thermal Weibel level 

So far we have been dealing with linear growth of the Weibel 
instability. In order to obtain its saturation level one needs to 
investigate the nonlinear evolution of the Weibel instability. 
Here we ask for how long it takes the Weibel instability to 
grow from thermal fluctuation level until reaching any mea- 
surable magnetic field strength. 

To estimate the thermal Weibel level we refer to thermal 



fluctuation theory of non-magnetic plasmas (Sitenko 1967 



[19821 |Akhiezer et al\ [19751 |Yoonl |2007b| l. (Working in the 
Weibel low frequency limit implies that in their general ki- 
netic expressions for the field fluctuations the variable z ~ w 
must ultimately be set to zero.) Recently the Weibel ther- 
mal level in an isotropic thermal has been calculated explic- 



itly (Yoon 2007b i. There a very simple expression for the 
wave number dependence of the spectral energy density of 
the magnetic field fluctuations has been obtained as 



(|b(0,k)| 2 ) th cxl/fc 3 Ai„ 



\d = Debye length 



(7) 



According to this expression the magnetic fluctuations 
maximise at long wave lengths k — > 0. In order to avoid 
and understand this effect, we assume that the thermal back- 
ground plasma in the current sheet has a weak thermal 
anisotropy T± ^ Tu. The direction of this anisotropy is ar- 
bitrary. Hence, for simplicity we assume that the 1 1 direction 
refers to the direction of the beam (perpendicular to the cur- 
rent sheet), while _L refers to the direction of current flow. 
This choice is not unreasonable in view of the possibility that 
the mere presence of the current may correspond to a thermal 
anisotropy. In equilibrium the background plasma then obeys 
a bi-Maxwellian velocity distribution 



/o(«||,«l) : 



(m e /27r)i 

t ± Jt\\ 



-exp 



2Tl 



2T» 



(8) 



In the classical limit of thermal fluctuation theory the mag- 
netic spectral energy density (bibj)koj in an isothermal 
plasma is determined from Eq. (2.52) in Sitenko ( 1967} as 

IhiTir ( < A-jfeA Imh U\k.S 



(bibj) 



k 2 



\e ± -n 



212 



(9) 



In a thermally anisotropic plasma one must replace the tem- 
perature T by the effective temperature T ± T\\/(T ± +Ty) in 
this expression. Moreover, the transverse dielectric function 
ej_(w,k) in the anisotropic case is given by 

, ,2 



c x (w,k) = l-- 



1- 



T 



l-$(z)+i7T2ze 



(10) 



Here z = u)/"\/2kv e u is a variable that vanishes with u) — > 0, 
n 2 = (kc/u)) 2 = e± is the refraction index of transverse fluc- 
tuations, and the real function $(2) ~ 2z 2 for z <C 1. 
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With the help of these expressions the spectral energy den- 
sity (Ibl 2 )/^ can be brought into the form 



y/^(cA; e ||)?||[fl/(l + l/g)] K e- z 
u e [K 2 -w 2 + l-0(l-$)] 2 +7rl9 2 z 2 e- 22; 



(11) 



where we defined k = k\ e ,£j = uj/uj e ,9 = T±/Tn, and A e = 
c/uj e is the background plasma electron inertial scale. This 
equation holds for small frequencies u ~ 0. Moreover, one 
now realises that for fixed ui = const and small k, using 
$(z) <~ 1 + 1 /2z 2 H — asymptotically for z>l, the spectral 
energy density becomes 



|b| 2 ) few oc 



KG 



[k 2 - Co 2 + 1 + 6/2z 2 } 2 + ir6 2 z 2 e- 2z2 



(12) 



which, as required, vanishes exponentially in the limit k <~ 
k — > 0, which holds for all frequencies to. 

The spectral energy density of the Weibel field is obtained 
in the limit z = u = 0, which yields 



|b| 2 ) fe0 : 



2 M0 yV2( C K||)T||[g/(l + l/g)]« 
w e " " [k 2 + 1-6»] 2 



(13) 



The isotropic case 9 — 1 reproduces the k~ 3 scaling with 
wave number which diverges at k — > 0. Retaining the temper- 
ature anisotropy the spectral energy density of Weibel fluctu- 
ations vanishes both for k — > and k — > oo. It, hence, has a 
maximum for 6* < 1 (i.e. Xj_ < Ty) at 



The value of this maximum is 



\ h \ 2 )k m a - 



3^ ^/2(c/v e]l )me/(i + i/e)} 



Suj e (1-9) 
which, for small anisotropies simplifies to become 



(14) 



(15) 



(16) 



The condition T± < T\\ corresponds to Weibel equilibrium. 
In the opposite case Tj_ > Tj| the thermal fluctuations ex- 
plode at wavenumber 



K 6 ■ 



(17) 



which indicates that in the vicinity of Kg the plasma is unsta- 
ble with respect to the ordinary anisotropic Weibel instabil- 
ity. The unstable fluctuation wavenumber Kg becomes small 
and the fluctuations long-wavelength for small anisotropies 
6 > 1. Then the plasma spontaneously pumps energy into 
the zero-frequency magnetic fluctuations while leaving the 
thermal equilibrium state. This would be the case when the 
sheet current generates sufficient thermal anisotropy to spon- 
taneously produce magnetic field fluctuations in the direc- 
tion perpendicular to the current; once this happens, a self- 
consistent magnetic field component will be created in the 



z-direction which could cause set-on of reconnection. Here, 
this particular interesting case will not be further considered. 

In a current sheet of density TV <~ 10 6 m~ 3 , tempera- 
ture T~ 0.1 keV which implies electron inertial and De- 
bye lengths A e w 6 km and \d ss 7.5 m respectively, pre- 
sumably corresponding to conditions in the near-Earth mag- 
netotail current sheet, the maximum Weibel spectral energy 
density becomes 



|b| 2 ) fcm0 «3x!0 



-26 J cV 



V 2 s- 



2„3 



-3 m 



(18) 



where T c y is the background plasma electron temperature in 
eV, and N cm s is the plasma density in cm~ 3 . 



4 Growth time in the magnetotail current layer 

We are interested in the time required for the Weibel instabil- 
ity to grow to measurable guide magnetic field values under 
conditions in the magnetotail reconnection region. Assuming 
an electron 'beam' density w 10 5 m -3 corresponding to 
Nb/N ~ 0.1, beam velocity V& <~ 30 kms -1 , angular beam 
size Ac* < 10°, the maximum growth rate of the Weibel in- 
stability in the tail current sheet becomes 



7W,max < 2 S 1 



(19) 



The Weibel instability picks up the thermal spectral energy 
density at its wave number of maximum growth fc w max > 
[1 + 2N/N b ]/X 2 eb which it amplifies. With the above num- 
bers the Weibel wavenumber in the magnetotail current sheet 
becomes 



k 



W.max^eb 



A eb >\/5 = 2.24 



(20) 



which corresponds to a maximum wavelength of Aw, max < 
2.8A e h in the magnetotail. 

The thermal fluctuation energy density at maximum grow- 
ing wave number is 



|b| 2 )fc w o< 



2^ yfiip(c/v e \\)T n [0/(l + l/6)]K W 



[«^ + l-*] s 



(21) 



where k w = fc w ,maxA e ~ (\ e /\ eb )[l + 2N/N b ] 2 > 0. In the 
isotropic case this simplifies to 



|b|W ^-4-1-4x10-- ^ 

W e V 2 V e K w HI 



(22) 



Because the maximum Weibel wavenumber kw 7^ is finite, 
this expression is free of divergence. 

We can now make use of the quasilinear instability growth 
prescription for the time evolution of the spectral energy den- 
sity 



(|b(i,k,0)| 2 )«(|b 2 (k,0)|) wth exp(2 7w i) 



(23) 
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where on the left is the linearly growing spectral density at 
time t, and on the right the thermal level from where the in- 
stability starts growing. No nonlinear saturation effects are 
included at this time. An observable magnetic field in the 
magnetotail should be roughly of order b <~ 1 nT, a value 
which one may use to find the required growth time tw from 



1 



|b 2 (k,T W )|) 
2 7w '"(|b 2 (k,0)|) w> th 



In 



(24) 



Inserting for the thermal energy density (|b 2 (k,0)|) w th 
from (|b| 2 )fe w0 and using the spectral energy density of a 
b ~ 1 nT magnetic field, {b{ nT ) fc0 ~ 4.3 x lCT 12 V 2 s 3 /m, the 
typical growth time to reach this magnetic field level be- 
comes 



TW,: 



(10-20) S 



(25) 



corresponding to 20-40 e-folding times. This time is in the 
range of typical onset times of substorms and thus not in con- 
tradiction to the onset of reconnection in the magnetotail. 

Since the maximum growing Weibel wave number is not 
precisely known (only a lower limit has been determined 
above), we can also use the maximum thermal fluctuation 
level at zero frequency which we have calculated for the 
anisotropic case. With its help we find for the thermal mag- 
netic pressure 



(Ibl 2 ) 



W,th 



2^o 



1.7 x 10" 27 (l-6») JriT 3 



0<|b|% 



(26) 



where we have used that, as has been argued above, the 
Weibel instability wave number k\ e > 1 is by order of mag- 
nitude comparable to the electron inertial scale A e . This ther- 
mal pressure corresponds to a thermal Weibel magnetic field 
strength 



(b 2 )w,th| 



5 ~6.5x lO^Vl-fl T 



(27) 



It is approximately this magnetic fluctuation level from 
where the Weibel instability in the current sheet starts to 
grow. Again, the quasilinear evolution of the spectral density 
can be used. However, maximum growth is practically inde- 
pendent of wave number k, and the same expression holds 
as well for the square of the maximum unstable Weibel- 
magnetic field. If we assume that the instability can freely 
grow until reaching observable values before nonlinear ef- 
fects come into play, the typical growth time for the magnetic 
field becoming observable is 



1 



T W ! 



2j 



In 



w 



|b(r w )| 2 ) 
Ilbl 2 ) 



W.th 



(28) 



We again demand that |b(r\y)| ~ 1 nT in order to be ob- 
servable in the magnetospheric tail plasma sheet. Since the 
calculation is anisotropic we assume the presence of a week 



anisotropy 1 — 6 <~ 0.01 only to again find a typical growth 
time of 



T W (lnT)wl0s 



(29) 



which corresponds to ^20 e-folding times. This is about 
the same time as that estimated above. Hence, the two ap- 
proaches are mutually consistent. 

During this growth time the beam traverses a distance of 
roughly ~ (600-1000) km 3> A e . This implies that during 
the growth time of the Weibel field the electron beams cross 
the entire electron inertial region, i.e. the inner field-free re- 
gion of current flow, thereby a posteriori justifying our as- 
sumption of symmetric counterstreaming electron beams in 
the magnetotail. 



5 Discussion and conclusions 

The result of this investigation is that inside the magnetic 
field-free electron inertial region (of transverse size of a few 
A e ) in the centres of thin - possibly reconnecting - current 
sheets, the inflow of electrons into the sheet from its two 
sides (in the geomagnetic tail from the lobes) may well be ca- 
pable of self-consistently generating a weak magnetic guide 
field via the non-magnetic Weibel instability. Given suffi- 
cient time, such weak guide fields will evolve in the very 
centre of the current sheet in the electron inertial region only. 
(Similar weak quadrupolar guide fields are caused in the 
much more extended ion inertial zone by the Hall effect, but 
these fields are found outside the current sheet centre, do not 
penetrate into the electron inertial zone and, thus, are only 
indirectly involved into the reconnection process.) 

The condition for the Weibel mechanism to work is that 
the ambient current sheet electrons have become nonmag- 
netic earlier than inert. This is the case that is realised in the 
centre of the magnetotail current sheet where the electron gy- 
rodradius exceeds the electron inertial length. In terms of the 
electron thermal speed this means that initially 



v e > yJmi/m e V A o 



(30) 



where Vao is the Alfven speed outside the current layer based 
on the external field. For an initial Harris current layer one 
writes this as a condition on the vertical coordinate z 



— < sinh 1 
A 



V A0 



(31) 



In the magnetospheric tail we have Vao ~ 10 6 ms _1 , and in 
the plasma sheet T e <~ 10 2 eV which yields that z < 0.1 A. 
In the inner region of the tail current this condition is clearly 
satisfied before entering the electron inertial region because 
A e - 0.02A,, and A - few A,. 

The Weibel-guide field B g is limited to be weaker than the 
ambient external magnetic field B - In the symmetric mag- 
netospheric tail current sheet we find that it may reach up to 
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< 10% of the ambient field in strength in a reasonably short 
linear e-folding time unless nonlinear effects set on earlier to 
saturate the field on a much lower level. Such nonlinear pro- 
cesses have not been investigated here. They can in principle 
be treated by assuming a stationary final state and calculating 
the thermal saturation level of the beam driven Weibel insta- 
bility. This, however, makes sense only if it is assumed that 
the stationary state will be reached in a time shorter than the 
onset of reconnection. The two processes cannot be consid- 
ered separately, since reconnection heavily depends on the 
presence of the (here self-consistently generated) guide field 
and will necessarily be fast, as has been shown by various 
numerical simulations that include the presence of even weak 
guide fields (see e.g., |Pritche"tt||2005[|Cassak et al.||2007) . 

Guide fields are important in the dynamics of the current 
sheet and in particular for reconnection. They have been 
observed in the geomagnetic tail in various cases (see e.g., 
Nakamura et al. 20081). They partially re-magnetise the elec- 



tron plasma in the central current region. Pointing along the 
electric field that drives the current, they cause particle accel- 
eration, which amplifies the current, generates energetic par- 
ticles, and generates a number of secondary effects that affect 
the stability of the plasma in the current sheet. The role of 
guide fields in collisionless magnetic reconnection and the 
various effects it may cause have in the past decade been 
thoroughly investigated in numerical simulations both for 
guide fields perpendicular (e.g., IRicci et al.[|2003 
and parallel (e.g 



Pritchett 2005 



2004a b) 



2006; Cassaketal. 2007 



Egedal et al. 2008 Le et al. 2009 ) to the sheet current. We 



may therefore refer to the published literature where the ex- 
istence of guide fields has been imposed. 

The interesting point we make here is that the Weibel 
instability provides a selfconsistent mechanism for produc- 
ing such guide fields in the field free current region. Since 
the plasma there is collisionless on quite a long time scale, 
the existence of guide fields is a nontrivial problem. Exter- 
nal fields can penetrate the plasma only on very long dif- 
fusion time scales, much longer than permitted by the ob- 
servation of reconnection onset. Hence, a mechanism that 
generates internal guide fields is highly welcome. The non- 
relativistic Weibel instability may provide such a mechanism. 
The Weibel fields are comparably small-scale forming mag- 
netic vortices of alternating polarity. On the scale of their 
wavelength, they cause a modulation of the magnetic field 
and consequently a spatial modulation of the current sheet. In 
the general case when the current sheet is not homogeneous, 
the guide fields may cause a three-dimensional structuring of 
the current sheet. 

The calculations presented here applied to the case when 
the conditions are homogeneous along the current, such that 
b • z = 0. In the general case, when homogeneity is not given, 
the Weibel magnetic field may possibly develop a non-zero 
field component b z = i„wz in the direction normal to the 
current sheet. Such a local normal magnetic field component 
is equivalent to a seed-X-point field and may ignite recon- 



nection in a similar way as assumed for initial condition in 
numerical simulations. 

Proof of this conjecture requires a full three-dimensional 
investigation of the non-magnetic Weibel instability under 
conditions when the current sheet is locally compressed (as 
sketched in Figure |2| in both x and y directions. If this 
conjecture turned out to be true, the local collisionless gen- 
eration of Weibel magnetic fields in thin current sheets in- 
side the electron inertial region \z\ < few A e would pro- 
vide a natural mechanism for initiating reconnection. The 
growth time of reconnection in such a case is the sum of 
the growth time tw of the Weibel instability starting from 
thermal level (b 2 ) t h until the normal Weibel field compo- 
nent 6„w becomes strong enough to modulate the ambient 
magnetic field, plus the well-explored growth time r rec of 
reconnection starting from an initial magnetic seed configu- 
ration. The latter is accessible through (and from existing) 
numerical simulations. 

Finally, we remark that under anisotropic conditions with 
Tj_ > Tji the Weibel thermal level diverges, which indicates 
instability for a certain wavenumber kq that depends on 
anisotropy 9 > 1. Such conditions can, in principle, cause 
spontaneous growth of a magnetic field component b z per- 
pendicular to the current sheet which would directly act as a 
seed field for spontaneous reconnection. 

Acknowledgements. This research is part of a Visiting Scientist 
Programme at ISSI, Bern. Hospitality of the ISSI staff and directors 
is thankfully recognised. 



References 

Achterberg A and Wiersma J: The Weibel instability in relativis- 
tic plasmas. I. Linear theory, II. Nonlinear theory and stabiliza- 
tion mechanism, Astron. Astrophys. 475, 1, doi: 10.1051/0004- 
6361:20065365,2007. 

Akhiezer A I, Akhiezer I A, Sitenko R V, and Stepanov K N: Non- 
Linear Theory and Fluctuations, in: Plasma Electrodynamics, 
Vol. 2, pp. 1 16-142, Pergamon Press, Oxford, 1975. 

Baumjohann W, Roux A, Le Contel O, Nakamura R, Birn J, 
Hoshino M, Lui A T Y, Owen C J, Sauvaud I-A, Vaivads A, 
Fontaine D, Rounov A V: Dynamics of thin current sheets: Clus- 
ter observations, Ann. Geophys., 25, 1365-1389, 2007 

Califano F, Cecchi T and Chiuderi C: Nonlinear kinetic regime of 
the Weibel instability in an electron-ion plasma, Phys. Plasmas 
9, 451, doi: 10.1063/1.1435001, 2002. 

Cassak P A, Drake J F and Shay M A: Catastrophic onset of fast 
reconnection with a guide field, Phys. Plasmas 14, 054502, doi: 
10.10 63/1.2734948, 2007. 

Drake J F, Swisdak M, Cattell C„ Shay M A, Rogers B N and Zeiler 
A: Formation of electron holes and particle energization dur- 
ing magnetic reconnection, Science 299, 873, doi: 10. 1126/sci- 
ence. 1080333, 2003. 

Drake I F, Shay M A, Thongthai W and Swisdak M: Production 
of energetic electrons during magnetic reconnection, Phys. Rev. 
Lett. 94, 095001, doi: 10.1103/PhysRev Lett.94.095001, 2005. 



R.A. Treumann, R. Nakamura, and W.Baumjohann: Weibel fields in thin current layers 



9 



Drake J F, Swisdak M, Che, H and Shay, M A: Electron acceleration 
from contracting magnetic islands during reconnection, Nature 
443, 553, doi: 10.1038/nature 05116, 2006. 

Egedal J, Fox W, Katz N, Porkolab M, 0ieroset M, Lin R P, 
Daughton W, and Drake J F: Evidence and theory for trapped 
electrons in guide field magnetotail reconnection, J. Geophys. 
Res. 113, A12207, doi: 10.10 29/2008JA013520, 2008. 

Ezawa Z F: Quantum Hall Effects, World Scientific, Singapore, 
2000. 

Fonseca R A, Silva L O, Tonge T W, Mori W B and Dawson J M: 
Three-dimensional Weibel instability in astrophysical scenarios, 
Phys. Plasmas 10, 1979, doi: 10.10 63/1.1556605, 2003. 

Fujimoto M, Nakamura M S, Shinohara I, Nagai T, Mukai T, Saito 
Y, Yamamoto T and Kokubun S: Observations of earthward 
streaming electrons at the trailing boundary of a plasmoid, Geo- 
phys. Res. Lett. 24, 2893, doi: 10.1029/97GL02821, 1997. 

Haaland S E, Sonnerup BUO, Dunlop M, Balogh A, Georgescu E, 
Hasegawa H, Klecker B, Paschmann G, Puhl-Quinn P, Reme H, 
Vaith H and Vaivads A: Four-spacecraft determination of mag- 
netopause orientation, motion and thickness: comparison with 
results from single-spacecraft methods, Ann. Geophys. 22, 1347- 
1365, 2004. 

Haaland S E, Sonnerup BUO, Dunlop M, Georgescu E, Paschmann 
G, Klecker B and Vaivads A: Orientation and motion of a dis- 
continuity from Cluster curlometer capability: Minimum vari- 
ance of current density, Geophys. Res. Lett. 31, L10804, doi: 
10.1029/2004GL020001, 2004. 

Le A, Egedal J, Daughton W, Fox W and Katz N: Equations of State 
for Collisionless Guide-Field Reconnection, Phys. Rev. Lett. 
102, 085001, doi: 10. 1103/PhysRevLett.l02.085001, 2009. 

Medvedev M V and Loeb A: Generation of magnetic field in the 
relativistic shock of gamma-ray burst sources, Astrophys. J. 526, 
697, doi: 10.1086/308038, 1999. 

Nagai T, Shinohara I, Fujimoto M, Hoshino M, Saito Y, Machida 
S and Mukai T: Geotail observations of the Hall current system: 
Evidence of magnetic reconnection in the magnetotail, J. Geo- 
phys. Res. 106, 25929, doi: 10.1029/2001JA900038, 2001. 

Nakamura R, Baumjohann W, Runov A, Volwerk M, Zhang T L, 
Klecker B, Bogdanova Y, Roux A, Balogh A, Reme H, Sauvaud 
J-A and Frey H U: Fast flow during current sheet thinning, Geo- 
phys. Res. Lett. 29, 2140, doi:10.1029/2002GL016200, 2002. 

Nakamura R, Baumjohann W, Asano Y, Runov A, Balogh A, 
Owen C J, Fazakerley A N, Fujimoto M, Klecker B and 
Reme H: Dynamics of thin current sheets associated with 
magnetotail reconnection, J. Geophys. Res. Ill, A11206, 
doi:10.1029/2006JA011706 A11206, 2006. 

Nakamura R, Baumjohann W, Fujimoto M, Asano Y, Runov 
A, Owen C J, Fazakerley A N, Klecker B, Reme H, Lucek 
E A, Andre M and Khotyaintsev Y: Cluster observations 
of an ion-scale current sheet in the magnetotail under the 
presence of a guide field, J. Geophys. Res. 113, A07S16, 
doi: 10. 1029/2007JA0 12760, 2008. 

Oieroset M, Phan T D, Fujimoto M, Lin R P and Lepping R P: In 
situ detection of collisionless reconnection in the Earth's magne- 
totail, Nature 412, 414, doi: 10.1038/35086520, 2001. 

Paschmann G, Haaland S, Sonnerup BUO, Hasegawa H, 
Georgescu E, Klecker B, Phan T D, Reme H and Vaivads A: 
Characteristics of the near-tail dawn magnetopause and bound- 
ary layer, Ann. Geophys. 23, 1481-1497, 2005. 



Pritchett P L: Onset and saturation of guide-field magnetic recon- 
nection, Phys. Plasmas 12, 062301, doi: 10. 1063/1.1914309, 
2005. 

Pritchett P L: Relativistic electron production during guide field 
magnetic reconnection, J. Geophys. Res. Ill, A10212, doi: 
10.1029/2006JA011793, 2006. 

Ramos J J, Porcelli F and Verastehui R: Driven reconnection about 
a magnetic X-line with strong guide component, Phys. Rev. Lett. 
89, 055002, doi: 10.1103/PhysRevLett.89. 055002, 2002. 

Retino A, Sundkvist D, Vaivads A, Mozer F, Andre M and Owen C 
J: In situ evidence of magnetic reconnection in turbulent plasma, 
Nature Phys. 3, 236, doi: 10. 1038/ nphys574, 2007. 

Ricci P, Lapenta G and Brackbill J U: Electron acceleration and 
heating in collisionless magnetic reconnection of a guide field, 
Phys. Piasmas 10, 3554, doi: 10.1063/1.1598207, 2003. 

Ricci P, Brackbill J U, Daughton W and Lapenta G: Collisionless 
magnetic reconnection in the presence of a guide field, Phys. 
Plasmas 11, 4102, doi: 10.1063/1.1768552, 2004a. 

Ricci P, Brackbill J U, Daughton W and Lapenta G: Phys. Plasmas 
11, 4489, doi: 10.1063/1.1778744, 2004b. 

Runov A, Nakamura R, Baumjohann W, Treumann R A, Zhang T 
L, Volwerk M, Voros Z, Balogh A, Glassmeier K-H, Klecker 
B, Reme H and Kistler L: Current sheet structure near magnetic 
X-line observed by Cluster, Geophys. Res. Lett. 30, 1579, doi: 
10.1029/2002GL016730, 2003. 

Silva L O, Fonseca R A, Tonge J W, Mori W B, Warren and Daw- 
son J M: On the role of the purely transverse Weibel instabil- 
ity in fast ignitor scenarios, Phys. Plasmas 9, 2458, doi: 10. 
1063/1.1476004,2002. 

Sitenko A G: Electromagnetic Fluctuations in Plasma, Academic 
Press, New York, 1967. 

Sitenko A G: Fluctuations & Non-linear Wave Interactions in Plas- 
mas, International Series in Natural Philosophy Vol. 107, Perga- 
mon Press, Oxford, 1982. 

Sonnerup BUO: Magnetic Field Reconnection, in: Solar System 
Plasma Physics, Vol III, pp. 45-108, eds. L T Lanzerotti, C F 
Kennel and E N Parker, North-Holland, New York, 1979. 

Sundqvist D, Retino A, Vaivads A and Bale S D: Dissipation in tur- 
bulent plasma due to reconnection in thin current sheets, Phys. 
Rev. Lett. 99, 025004, doi: 10.1 103/PhysRevLett.99.025004, 
2007. 

Vaivads A, Khotyaintsev Y, Andre M, Retino A, Buchert C, Rogers 

B N, Decreau R, Paschmann G and Phan T D, Phys. Rev. Lett. 

93, 105001, doi: 10. 1 103/PhysRevLett.93. 105001, 2004. 
Weibel, E S: Spontaneously growing transverse waves in a plasma 

due to an anisotropic velocity distribution, Phys. Rev. Lett. 2, 83, 

doi: 10.1103/PhysRevLett.2.83, 1959. 
Yoon P H: Relativistic Weibel instability, Phys. Plasmas 14, 02454, 

doi: 10.1063/1.2646285, 2007a. 
Yoon P H: Spontaneous thermal magnetic field fluctuations, Phys. 

Plasmas 14, 064504, doi: 10.1063/1.2741388, 2007b. 
Yoon P H and Davidson R C: Exact analytical model of the classical 

Weibel instability in a relativistic anisotropic plasma, Phys. Rev. 

A 35, 2718, doi: 10.1 103/PhysRevA.35.2718, 1987. 



